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Abstract. The properties of discrete two-dimensional spin glasses depend strongly 
on the way the zero-temperature limit is taken. We discuss this phenomenon in the 
context of the Migdal-KadanofF renormalization group. We see, in particular, how 
these properties are connected with the presence of a cascade of fixed points in the 
renormalization group flow. Of particular interest are two unstable fixed points that 
correspond to two different spin-glass phases at zero temperature. We discuss how 
these phenomena are related with the presence of entropy fluctuations and temperature 
chaos, and universality in this model. 
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Since the beginning of the study of disordered systems with renormalization and 
scaling methods the question of universality and of the relevance of the realization of 
the disorder have been a key issue pQ. Why indeed should we accept that some abstract 
models are the archetype of a very broad class of systems if their behavior depends 
drastically on tiny details? The Edwards- Anderson [2] model is one of the models that 
are widely regarded as a prototype of disordered systems in statistical physics. It is an 
Ising model with disordered and competitive interactions and has been the source of 
many surprises and developments in the last thirty years [31 H] . 

Its two-dimensional (2D) version, one of its simplest settings, has very special 
properties with respect to universality. It is now agreed that the spin-glass phase 
exists only at zero temperature EJ [7], where the spin-glass susceptibility diverges. 
However, the behavior of the model seems to depend drastically on microscopic details 
and in particular discrete and continuous couplings leads to different properties [El |9] 
at zero temperature. On the other hand, for nonzero temperatures, strong evidence 
for universal critical behavior has been observed [TQl [HI [12]. This rises questions on 
the very nature of universality, if any, in strongly disordered bidimensional systems. 
Why such a difference? What are the mechanisms behind this behavior? The key 
to understand these features lies in the difference between strictly zero and vanishing 
temperature [131 EH]- When temperature is finite, entropy fluctuations play a major 
role and, for large enough sizes, bring back discrete models to the continuous class 
[TO] . These mechanisms were further exploited in [15], with a special emphasis on the 
so-called temperature chaos effect [161 El 118] . 

In this Letter we study the behavior of continuous and discrete 2D spin glasses in 
the context of the Migdal-Kadanoff renormalization group (MKRG) [19] . We observe 
that indeed continuous and discrete 2D spin glasses are eventually associated with the 
very same physically relevant fixed point. However, there are also key differences that 
are the consequence of a cascade of two repulsive fixed points in the MKRG flow in 
the case of the discrete model. In a real-space picture, these phenomena are associated 
with two different temperature-dependent crossover length scales, so that two different 
zero-temperature spin-glass phases can be observed depending on the system size L and 
the temperature T. 

This Letter is organized as follows: first, we define the model and recall briefly 
the MKRG. We then study the MKRG flow and observe one attractive paramagnetic 
fixed point and two repulsive ones corresponding to the T = discrete and continuous 
classes, respectively. We discuss the crossover length scales between different regimes, 
the cascade of fixed points and its effects on effective critical exponents which may easily 
camouflage the fact that the critical behavior of the 2D spin glass is universal. 
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1. Two-dimensional spin glasses 

The Hamiltonian of the Edwards- Anderson Ising spin glass [2] is given by 

H = - 2J JijSiSj. (1) 

<i,j> 

The spins Si lie on a square lattice of size N = L 2 in two space dimensions and the 
interactions between the spins are between nearest neighbors. We shall consider two 
different models: in the discrete onaj, the interactions are chosen randomly as G {±1} 
while in the continuous, or Gaussian model, we shall use a Gaussian distribution of the 
Jij with zero mean and unit variance, instead. 

It is by now generally accepted (although there is no rigorous proof of this, but 
see[20]) that there is no spin-glass transition at finite temperature in both cases. Instead, 
the spin-glass susceptibility diverges when T — > and the spin-glass phase exists only 
at zero temperature. At any finite temperature, however, this implies the existence 
of an equilibrium length scale £ e q(^) beyond which the spin-glass correlation decays 
exponentially fast and the system is effectively in a paramagnetic state. For length 
scales below £ eq (T), the correlation function decays only as a power-law and the system 
has a spin-glass-like correlation. If one is looking at a system of size L <C £ eq (T), it thus 
looks very much as a spin glass (it has power-law correlations) and only for larger sizes, 
when L 3> Ceq(^) does the system finally look paramagnetic. 



2. Migdal-Kadanoff renormalization group 

Low-dimensional systems are often well described by the Migdal-Kadanoff renormaliza- 
tion group and 2D spin glasses are no exception [T9l I2TI |22| [23] . We shall thus work 
within the MKRG approach and study how different are 2D discrete and continuous 
spin glasses. For a disordered system such as a spin glass, the MKRG is a functional 
renormalization group method, as the quantity that is being renormalized is the distri- 
bution of couplings P(/3J). Starting from a given inverse temperature /3 = 1/T and an 
initial distribution of couplings P init (/3 J), MKRG allows to follow the (approximate) flow 
under renormalization (for a detailed description see for instance [2T1 1221 123]). There 
are many ways to write the MKRG recursion for 2D spin glasses. We shall follow [9] 
and use the recursion on a hierarchical lattice with 6 = 3 branches and s = 3 spins per 
branch yielding a model with effective dimension D = 2: 

r 9 

P G+ \PJ)= /n (dPiP G {PJi)) &W-H{PJi}i=i... 9 )\ (2) 
J 1=1 

with J 7 ({/5Ji} i=1 9 ) = atanh(tanh/3 J\ tanh/3 J 2 tanh/3 J 3 ) 

+ atanh(tanh / 9J 4 tanh / 9J5tanh/3J 6 ) 
+ atanh(tanh/3J7tanh/3Jgtanh/3Jg). 

X In this letter, discrete means discrete energy spectrum, thus not allowing for irrational discrete values 
for the coupling Jy as in [9]. 




Figure 1. Top: The domain-wall free energy Ai<bw = \J {(Jp — Ja.) 2 ) associated with 
an effective bond after G MKRG iterations (where the effective size of the system 
is L ~ 3 G ), in the low-temperature phase of the Gaussian (left) and binary (right) 
spin glasses. In the Gaussian case one is in the spin-glass phase where AFdw oc 3 , 
with 9 w —0.278 until the effective size of the system is larger than a crossover length 
£,eq(T) = T~ u where the system becomes paramagnetic and 9 = — oo. In the discrete 
case, the system is first in a phase where 9' — 0, until a first length scale £ c is crossed, 
and only at this point it decays with the continuous exponent 9 « —0.278. When 
the effective system is even larger, such that L 3> Lq oc T _I/cff , the system becomes 
again paramagnetic: there are therefore two different spin-glass phases in the zero- 
temperature limit. Bottom: the corresponding flow in the MKRG with the three 
different fixed points. Binary spin glasses are cascading through the three fixed points. 



This recursion can be implemented trivially using population dynamics. It is also 
sometimes convenient to see the MKRG as an exact solution on a hierarchical lattice [24J . 
In this case one interprets the "time" in the RG flow (that is, the index G in Eq. [2]) as 
related to the effective size of the system via L = 3 G . This allows to characterize length 
scales in the system, and it also provides an exact realization of the droplet /scaling 
theory pSJ [T71 [18] , which is, at least for the 2D model, accepted to be a good description. 
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MK 


2D 


e 


-0.278 


-0.287 [EH] 


v = -1/9 


3.597 


3.48 


d s 


1 


1.28 [261127] 


C=f -9 


0.778 


0.927 



Table 1. Gaussian couplings: critical exponents from the MKRG approximation and 
the 2D square lattice. 

3. Gaussian couplings 

When starting from a Gaussian distribution Pi n it(/3J) at low temperature, the MKRG 
flows to a symmetric slightly non-Gaussian distribution with zero mean and a variance 
decaying with G such that (J 2 ) oc L d , where 9 ~ —0.278. At strictly T = 0, this is 
the strong-coupling fixed point describing the continuous spin-glass phase, but at any 
finite temperature T the fact that 9 < implies that the bonds are not robust and that 
the spin-glass phase does not survive (and that the fixed point is unstable). This is 
precisely what happens once enough recursions are performed and the MKRG flows to 
a trivial delta function exponentially fast: this is the stable paramagnetic pseudo- fixed 
point where formally 9 = — oo (see Fig. [Hand Fig. [2j) 

Following the standard droplet interpretation, scaling relations can be obtained 
by realizing that the thermodynamic behavior is dominated by large excitations (the 
droplets) of size £ and 0(l)-energy that can be found with a probability scaling as £ e . 
The probability that two spins at distance r are correlated thus decays as r 9 and this will 
be on the order of the thermal fluctuations when r reaches the equilibrium length scale 
£ eq oc T l l e . In other words the correlation function behaves as CG auss (r, T) oc r e e~ T ^ c ^ T " > 
which indicates that the critical exponent for the correlation length is simply v — —1/9. 
Indeed, when performing the MKRG at finite T, we observe that when the effective size 
L is of the order of £ (which eventually happens for any positive temperature, provided 
we do enough recursions), the value of 9 drops from —0.278 to — oo. This is the sign 
that we are now in the paramagnetic phase and that there are no correlations beyond 
the scale £ e q(T) (see again Fig. CQ. 

We recall all the MKRG exponents for the Gaussian model, together with the 
estimated ones from the 2D square lattice in Table [TJ The values are very close, showing 
that the critical properties are well captured by the MKRG approximation. 

4. Discrete spin glasses 

The discrete ± J model displays a more puzzling phenomenology. It was first realized 
that strictly at zero temperature, one finds 9' = (we shall use a prime for all exponents 
in the discrete model that are related with this fixed point, that we will refer to as the 
discrete fixed point), both in MKRG [9] and in the 2D model [8], suggesting a different 
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MK 


2D 


9' 








v' = 1/C 


2.597 


1.43 ...1.83 




0.77 


1.095 ...1.395 [28] 


/*/ _ d s 

S 2 


0.385 


0.5 ...0.7 


^cff = f + V> 


6.195 


4.9 ...5.5 



Table 2. Discrete ±J couplings: critical exponents associated with the "discrete" 
fixed point in the MKRG approximation and the estimated ones from the 2D square 
lattice. The exponent v e s describes the divergence of the spin-glass correlation length 
using a scaling in the system size L and the bare temperature T . 

universality class with respect to Gaussian disorder. However, it was soon realized 
that a small perturbation (such as discrete, but fractional couplings) were enough to 
destabilize this fixed point, and the MKRG was then flowing to the continuous one [9] 
where 9 = —0.287. Subsequent Monte-Carlo simulations of the 2D model indicated that 
some observables were in the same universality class in both the discrete and continuous 
mo dels [TO]. This rises questions: how one goes from different results at zero temperature 
in the discrete and the continuous model to the same ones at finite T? We will see that 
there are actually two different unstable fixed points, corresponding to two different 
spin-glass phases, that can be observed depending on how the T = limit is taken, and 
a stable pseudo-fixed point corresponding to the paramagnetic phase. 

4-1. The zero-temperature discrete phase 

We first repeat the T = analysis of [9] and, indeed, find 9' = 0. That means that the 
P G (J) (which we rewrite in order to remove the dependence in /3 by taking formally the 
limit — > oo in Eq. [2]) converges to a nontrivial function that does not evolve under 
renormalization. This fixed point thus describes a zero-temperature spin-glass phase, 
whose properties are different from the ones seen in the Gaussian model. 

The nontrivial phenomenon, however, is that the flow is not going directly to the 
paramagnetic pseudo-fixed point, but instead it is first governed by the continuous fixed 
point, and then finally flows to the paramagnetic one (see Fig. [TJ implying a nontrivial 
cascade in the MKRG flow. 

4-2. Entropic fluctuations and temperature chaos 

The reasons behind the instability of the discrete fixed point are, as first discussed by 
[10] . entropic fluctuations that at any finite temperature make the zero-temperature 
and the vanishing-temperature limits different (a very similar phenomenon appears in 
diluted 3D spin glasses [U], where the critical dilution calculated by first setting the 
temperature to zero and then taking the system size to infinity is strictly larger than 
when first taking the system size to infinity at finite temperature and then taking the 
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Figure 2. (Color online) Effective 9 and fractal dimension of domain walls d s for the 
Gaussian and ± J models as a function of the number of generations G in the MKRG 
for different temperatures T. The length scale probed by the MKRG is L = 3 G . For 
the discrete model, the rescaling using the exponent v' = allows to superimpose 
the curves in the crossover region from the discrete to the continuous fixed point, while 
using i/ c ff = v + v 1 one can superimpose the crossover from the continuous fixed point 
to the paramagnetic region. 



temperature to zero). This is in fact deeply connected with the notion of temperature 

chaosnni m mi nil m us una. 

Equilibrium states of spin glasses are sensitive even to very small perturbations. 
Let us repeat here the classical thermodynamic argument of [161 HE]- Consider a 
spin glass in the scaling/droplet approach: take two equilibrium states at temperature 
T < T c differing by a very large droplet of characteristic size t. Then the two states 
have free energies that differ by AF(T) = AE — TAS ~ Y(T)£ 9 , where T(T) is 
the energy stiffness coefficient. When one changes the temperature by AT, then 
AF(T + AT) ~ AE — (T + AT)AS so that AF(T + AT) w T(T)^ - AT AS. In 
this phenomenological approach, the entropy difference is associated with the droplet 's 
surface so that AS" has a random sign and a typical mag nitude a(T)£ d ^ 2 , where <x(T) 
is called the entropy stiffness and d s is the fractal dimension of the droplet's surface. If 
d s /2 > 6, which follows from droplet theory, then AT(T + AT) can change sign between 
T and T + AT for length scales greater than 

with ( = d s /2 — 9. This can be checked in the MKRG using the correlation between the 
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Figure 3. (Color online): Temperature chaos in 2D spin glasses: we show ln(l — 
C(T, AT))/2 where C(T, AT) is the correlation between the couplings at temperature 
T and T + AT. In the Gaussian model, the system is characterized by a critical 
exponent £, while in the binary case one can observe either the exponent £ or 
depending on the values of L and T. 



effective couplings between the border spins at different temperatures 

_ (J(T)J(T + AT)) 
C[I ' AI} - (J(T))(J(T + AT)Y (4) 

Besides the scaling in L/£ c , it is also expected|30] that for small values of AT one 
has 1 — C oc L 2C - . As shown in Fig. [3] we find that indeed one can rescale the curves and 
observe the predicted slope using ( = 0.78 for the Gaussian system. Using the definition 
of C, this shows that d s — 1 in MKRG, as is widely accepted. For the discrete case, 
however, we first find a (' ~ 0.385 in the slope instead. This was to be expected since 
9' ^ 9. In turns, this indicates that d' s ~ 0.77. This is surprising, since in real systems d s 
cannot be lower than d — 1, and one has no choice but to accept it as a strange byproduct 
of the MKRG approximation. The point, however, is that once enough recursions are 
performed, the exponent in the slope is not given by the discrete fixed point exponent 
£' anymore, but by the continuous one (, and the points superimpose where rescaled 
by C : this is a clear sign that we are now effectively in the continuous spin-glass class 
(and this was in fact first noticed by [U]). We also find that d! s is smaller than d s in the 
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MKRG as well as on the square lattice model. 

In a real-space interpretation, one can picture the phenomenon by thinking of 
droplets "dressed" with entropic fluctuations: when the entropy plays a role, at finite 
temperature, the droplet's cost is given by an interplay between the values of J, which 
are discrete, but also by the entropic values proportional to T (given by some "spins 
fous" in effective zero field [13]) which is continuous: the free-energy cost of an excitation 
starts to be slightly continuous and the MKRG flow must leave the discrete fixed point 
to go to the continuous fixed point. The fact that this is due to entropy fluctuations, 
and therefore temperature chaos, is yet another illustration of the its importance in spin 
glasses. 

4-3. The continuous spin-glass regime 

For sizes larger than £ C (T), we thus expect the MKRG flow to be governed by continuous 
fixed point. Indeed, we recover the "continuous class" critical exponents d s and 9 once 
the number of iteration is large enough such that L > £ c : this can be seen directly in 
the behavior of the domain- wall exponent (Fig. [1] and Fig. [2]) where the three different 
regimes (discrete, continuous and paramagnetic) are observed successively. We have 
also checked, by scaling the values of 9 obtained after a given number of iterations G at 
temperature T, that the crossover from the discrete to the continuous fixed point arises 
at £ c , and is therefore controlled by the exponent of the discrete fixed point v' = 

4-4- The paramagnetic regime 

At any finite T, the flow must eventually end in the paramagnetic attractive pseudo- 
fixed point. Thus, according to the MKRG, there should be three different regimes: 
for sizes L < £ c we have 9' = 0. For £ c < L < L cq we have 9 = —0.278 while for 
L > L oq we are in the paramagnetic phase and 9 = — 00. The picture of the flow is 
one of a cascade of three fixed points (two unstable and a final stable one, see FigJT]). 
However, due to this cascade, the two repulsive fixed points combine to produce an 
apparent nontrivial exponent for the crossover finite-size length scale, which is given by 
^cfr = v + 1/C'> instead of v in the Gaussian problem. This can also be understood 
directly in the scaling picture: Let us consider the probability that two spins at distance 
r are coupled. Starting from size £ c = T 1 ^ it decays as p oc (r/£ c ) e . Indeed, we saw 
that the do main- wall free-energy is 0(1) for r = £ c . This indicates that this will be of 
order 0(T) for L oc T^-^T^ 1 ^' so that L eq oc T 1 ^^ 1 ^' . In other words the effective 
exponent using the system size and the bare temperature T as scaling variables is not 
z/, but z/ e fj = v + 1/C'- We will, however, see that the use of these scaling variables only 
camouflages the fact that the divergence of the spin-glass correlation length is simply 
given by the exponent v. 
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5. Universality 

Universality is a tricky problem in this model, precisely because of this nontrivial cascade 
of fixed points! Depending on the temperature T and the size L, one can be sensitive 
either to the "continuous" critical point, to the discrete one, or even a combination of 
the two in simulations, and it is important not to confuse the two. Once, however, the 
peculiar picture of the flow (Fig. [T|) is understood, it is clear that the critical fixed point 
governing the RG flow before the paramagnetic is reached is the same for both Gaussian 
and binary disorder. 

We shall now see how one can observe the continuous class in simulation of the 
discrete model. As argued in [3TJ [12] using bare (unrenormalized) couplings (e.g., 
temperature) as scaling variables is not a good idea in this model and leads to confusing 
results, but instead dressed (renormalized) couplings (e.g., Binder ratio) should be used. 
This was precisely the idea used in the simulations in 2D that observed strong evidence 
for universality [TOl [12]. A good examples of such a quantity is the Binder cumulant 



g(L,T) = |(3— j^yz ) of the overlap distribution P(q) and the spin-glass susceptibility 



X = N(q 2 ). These are classical parameters in simulations and they can be computed in 



MKRG as well (see [32] for details) which gives a good idea of what is to be expected 




in the 2D system on a square lattice at finite size. We display the scaling functions of 
the Binder ratio and the spin-glass susceptibility as introduced in [3T] for the Gaussian 
and the discrete disorder in Fig. [H These graphs allow one - amongst other things - to 
determine the critical exponents v and 77. The points obtained from the Gaussian and 
the discrete model superpose on a nontrivial universal curve, apart from the fact that the 
discrete model has an additional branch which is related to the discrete fixed point, but 
which, however, is irrelevant for the critical exponents observed in the thermodynamic 
limit. This is, indeed, precisely what is also observed in simulations of the 2D binary 
square lattice model [10]. It is interesting to note that, in contrast to the continuous 
fixed point, the relation 9 = —\jv does not hold for the discrete fixed point where 



Another universal quantity is the value of the Binder ratio at the critical point 
g c which for the Gaussian model is g c = 1. For the discrete model the situation is 
slightly more involved because of the cascade of fixed points in the RG flow. In the left 
panel of Fig. [5] we show the Binder ratio as a function of the system size (or number 
of generations G). It clearly shows the different scaling regimes. At small system sizes 
the discrete fixed point dominates the scaling behavior and increasing the system size 
initially leads to an increase of the value of the Binder ratio until, upon increasing the 
size even further, we enter the scaling regime of the continuous fixed point where the 
Binder ratio starts to decrease when the system is made larger. The maximal value 
that the Binder ratio reaches upon entering in the domain of the continuous fixed-point 
scaling depends on the temperature as can be seen on the right panel of Fig. [5] and close 
to zero temperature it scales linearly with a zero-temperature limit that is perfectly 
compatible with g c — 1, i.e., that the value of the Binder ratio at the critical point is 





0, but v' = 2.597. 



The nature of spin-glass phases in 2D discrete spin glasses 



11 



also the same in both models. The fact that g c = 1 and that 77 = also means that 
there is one pair of physically relevant ground states in both models. 




Figure 4. (Color online) In the left panel the scaling function of the Binder ratio 
g(L,T) is shown for that Gaussian and ±J model. The data collapse shows that they 
share same scaling function and hence the same critical exponent v. In the right panel 
the scaling function of the spin-glass susceptibility \ is shown. The data collapse on the 
physical branch of the scaling functions shows that the critical exponent r\ is universal. 



6. Discussion 



Using the MKRG approach, we have studied 2D spin glasses, in particular the model 
with discrete couplings and its two different spin-glass phases at zero temperature. In the 
thermodynamic limit, one can observe each of these phases depending on how the zero- 
temperature limit is taken: If one takes first the zero-temperature limit (or simply sets 
the temperature to zero right from the beginning as it is often done in numerical ground 
state studies) and then sends the size to infinity, the discrete fixed point is obtained. If 
however, one first sends the size L to infinity and then sends the temperature to zero, 
then entropy fluctuations (and temperature chaos) play a role and one is in fact observing 
a fixed point identical to the one reached with continuous couplings, strongly indicating 
that there is universality in 2D spin glasses. It is important to note that the physically 
relevant fixed point is the one which determines how quantities such as the spin-glass 
correlation length or the spin-glass susceptibility diverge as the temperature approaches 
the critical temperature T c = and this is in all cases the continuous-coupling fixed 
point. It is also important to note that physically the domain- wall free-energy exponent 
is 9 = \ jv = —0.278 also in the discrete coupling model if the thermodynamic limit is 
taken such that the system size is always kept larger than l c {T) which notably is not the 
case if one simply sets the temperature to zero in a finite-size scaling study (this leads 
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Figure 5. (Color online) In the left panel the Binder ratio g(G,T) is shown as a 
function of the number of generations G used in the construction of the hierarchical 
lattice. In the right panel the maximal value of the Binder ratio in shown as a function 
of temperature T. A linear extrapolation to zero temperature is very well compatible 
with universality, i.e., with g c = 1. 

to the well-known 9' = 0.) At finite size, the observed properties depend on the length 
scale under probe and the crossovers between the different regimes are associated with 
the chaotic length scale £ c and the paramagnetic one L eq . 

The study of the continuous unstable fixed point is made difficult in simulations 
because of the cascade in the MKRG flow. However, we have shown that when using 
proper RG invariant quantities, one can measure the exponents associated with the 
fixed points precisely. Let us finally point out that this situation, where the entropy 
changes drastically the properties at zero temperature, is not entirely new. In fact, this 
is precisely the effect that changes the location of the spin-glass transition in diluted 
3D spin glasses [14J from the naive percolation point to a nontrivial one: this shows how 
much entropic effects, and temperature chaos, matter in any renormalization or scaling 
scheme of finite-dimensional spin glasses [13] . In fact, such entropic effects also matter 
in mean- field spin glasses and optimization problems [33]. One has therefore always to 
be cautious in using ground-state properties if the entropy is not taken into account, 
because the presence of an additional unphysical fixed point might lead to wrong results. 

During the (long) preparation of this manuscript, we became aware of the work of 
|15j . who first stressed the role of temperature chaos. In fact the exponent 5.5 they find 
in 2D is precisely z/ e g in Table 2. The recent simulations of [12], as well as the older ones 
in [10], display strong indications supporting our analysis in 2D spin glasses on square 
lattices. 
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